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We introduce two pairs of stable cheapest nonconforming finite element space 
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nite element space; the pressure field is approximated by the piecewise constant 
function with mean-zero space eliminated. We show that two element pairs 
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tionship between them. Several numerical examples are shown to confirm the 
efficiency and reliability of the proposed methods. 
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1. Introduction 


In the simulation of incompressible, viscous fluid mechanics, the lowest- 
degree conforming element Pi x Pq or Qi x Pq produces numerically unsta¬ 
ble solutions in the approximation of the pressure variable [10]. In particular 
Boland and Nicolaides [3, 4] fully investigate for the pair Qi x Pq. The above 
simple pair does not satisfy the discrete inf-sup condition. Several successful 
finite elements satisfying this condition have been proposed and used. For in¬ 
stance conforming finite element spaces [2, 9, 25, 26] including the P 2 x Pq and 
P 2 X Pi (the Taylor-Hood element) elements [11, 13] and the MINI element [1]. 

Instead of conforming finite element spaces, the use of nonconforming fi¬ 
nite element spaces has been regarded as one of the simplest resolutions to 
the discrete inf-sup conditions: see [7] for simplicial elements with the Pi non- 
conforming element for the velocity approximation and the Pq element for the 
pressure approximation. For rectangular and quadrilateral elements, the use 
of nonconforming elements with four or five degrees of freedom with the pres¬ 
sure approximation by Pq element leads to stable element pairs for the Stokes 
equations [6, 8, 12, 14, 15, 18, 22, 24]. 

The use of Pi nonconforming quadrilateral element, whose local degrees of 
freedom are only 3, in the approximation of velocity fields with Pq approxi¬ 
mation to the pressure leads to unstable finite element spaces. An interesting 
question arises: what are the smallest rectangular/quadrilateral nonconform¬ 
ing element spaces to approximately solve the velocity fields combined with Pq 
approximation to the pressure? 

Recently, Nam et al. [20] introduced a cheapest rectangular element based 
on the Pi nonconforming quadrilateral element [21] by adding a globally one¬ 
dimensional bubble function space [8, 24] to the Pi x Pq pair on rectangular 
meshes. They show that the one-dimensional enhancement to the velocity space 
fulfills the discrete inf-sup condition whose constant depends on the mesh size 
h and provide several convincing numerical results with smooth forcing term. 
However, it has been questionable whether this one-dimensional modification 
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can lead to a stable cheapest element or not. 

The primary aim of this paper is to propose two stable cheapest finite el¬ 
ement pairs based on the Pi nonconforming quadrilateral element space and 
the piecewise constant element space. Our modification is still a globally one¬ 
dimensional enhancement to the velocity space enriched by adding a globally 
one-dimensional DS'5'T-type (or Rannacher-Turek type) bubble space based on 
macro interior edges. Equivalently we propose to modify the pressure space by 
eliminating a globally one-dimensional spurious mode with the velocity space 
unchanged from the Pi nonconforming quadrilateral element space (For a con¬ 
forming counterpart, see [10]). 

Indeed, these two finite element pairs are closely related. We show that 
the velocity solutions obtained by these two finite element pairs are identical 
while the pressure solutions differ only by a term ff[h) times the global discrete 
checker-board pattern. Thus, the stability and optimal convergence results for 
one finite element pair are equivalent to those for the other. 

It should be stressed that if the conforming bilinear element is used instead 
of our Pi nonconforming quadrilateral element with the same modification to 
the pressure space, the conforming bilinear element is still not stable (See Cor. 
5.1 and numerical results in Tables 4 and 5 in §5. 

Recently, the proposed elements are used to solve a driven cavity problem 
[17] and an interface problem governed by the Stokes, Darcy, and Brinkman 
equations [16]. 

The outline of this paper is organized as follows. In Section 2, the Stokes 
problem will be stated and the first finite element pair will be defined. In Section 
3, we define the second finite element pair and present a relationship between 
our two finite element pairs. Section 4 will be devoted to check the discrete inf- 
sup condition for our proposed finite element pairs by using a technique derived 
by Qin [23]. Finally, some numerical results are presented in Section 5. 
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2. The Stokes problem and the stabilization of pressure space 

In this section we will introduce a stable nonconforming finite element space 
pair for the incompressible Stokes problem in two dimensions. We begin by 
examining the pair of Pi nonconforming quadrilateral element and the piecewise 
constant element. Then a suitable minimal modification will be made so that 
uniform discrete inf-sup condition holds. 


2.1. Notation and preliminaries 

Let O C be a bounded domain with a polygonal boundary and consider 
the following stationary Stokes problem: 


— -t- Vp = 

f 

in fl. 

(2.1a) 

V • u = 

0 

in n, 

(2.1b) 

u = 

0 

on dil, 

(2.1c) 

where u = (iti, ^ 2 )^ represents the velocity 

vector, p the pressure, f = 

(/i,/2)^e 


H ^(fl) the body force, and > 0 the viscosity. Set 

LQ{fl) = {q G L^IQ.) I f g dx = 0}. 

Jn 

Here, and in what follows, we use the standard notations and definitions for the 
Sobolev spaces H'’(S'), and their associated inner products (•, Os.S; norms 11-115,5, 
and semi-norms | • |s,s. We will omit the subscripts s,S' if s = 0 and S' = H. 
Also for boundary dS of S, the inner product in L^{dS) is denoted by 
Then, the weak formulation of (2.1) is to seek a pair (u,p) G Hj)!!) x 
such that 

a(u,v) - &(v,p) = (f,v) VveHj(H), (2.2a) 

b{u,q) = 0 VgeLo(^)i (2.2b) 

where the bilinear forms a(-,-) : Hj(r2) x H(,(H) —R and &(-,-) ; Hj)!!) x 
Lq(H) —i- R are defined by 

a(u,v) = i/(Vu, Vv), 6(v,g) = (V-v,g). 
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Let ^ = {v G Hj(ri) I V-v = 0} denote the divergence-free subspace of Hj(0). 

Then the solution u of (2.2) lies in ^ and satisfies 

a(u,v) = (f, v) Vv G (2.3) 

2.2. Nonconforming finite element spaces 

In order to highlight our approach to design new finite element spaces, we 
shall restrict our attention to the case of = (0,1)^. Let {^y^h)o<h<i be a 
family of uniform triangulation of into disjoint squares Qjk of size h for 
j,k = 1, ■ • ■ ,N and fl = U^fe=i Qjk- denotes the set of all edges in ^h- 
Let Nq and iV® be the number of elements and interior vertices, respectively. 

Let Pj{Q) denote the space of polynomials of degree less than or equal to j on 
region Q. 

The approximate space for velocity fields is based on the Pi nonconforming 
quadrilateral element [5, 8, 21]. Set 

^nc,h _ j.^ g L^(fl) I v|q G Pi(Q) VQ G is continuous at the midpoint 
of each interior edge in 

and 

^ 1 , 0 ^ = {v G I V vanishes at the midpoint of each boundary edge in ^}. 

The pressure will be approximated by the space of piecewise constant functions 
with zero mean i.e., 

^0 = {« G I q\Q G PoiQ) VQ G ^4, dim = Nq - 1. 

It is known that the pair of spaces x cannot be used to solve the 

Stokes equations, as stated in the following theorem: 

Theorem 2.1 ([20]). Let {tPh)o<h<i be cl family of triangulations of LI into 
rectangles and set 

= {ph G 1 bh{vh,Ph) = 0 Vv/^ G 

where bh{vh,Ph) ■= ' "^h,Ph)Qi- Then dim('^'‘) = 1. Indeed, the ele¬ 

ments Ph G are of global cheeker-board pattern. 
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Denote by C-u a global checker-board pattern basis function with ||C/j|| = 1 
such that 


= Span{C/,}. (2.4) 

For simplicity, we assume that can be considered as the disjoint union of 
macro elements such that each macro element consists of 2 x 2 elements in 
For odd integers j and k, consider the macro element consisting of Qjki 

Qj,k+i^ Qj+i,k, and Qj+i^k+i, with {J,K) = {j,k). Denote by the macro 

triangulation composed of all such macro elements QjkS- Let G be 
the elementary checker-board pattern defined by 


-1 

1 

on Qfj, = 

Qj,k-\-l 

Q j-\-l,k-\-l 

1 

-1 


Qj,k 

Qj-\-l,k 


0 on D \ 








Q78 

Qs 8 

Qfi 


Ql{ 


QU 


Qn 

Qs 7 






Qj+i. 

:+l 


Qfl 


Qll 


Q^k 

Qj,k 

Qj-\-l.k 

Qli 










Qfi 


Qli 


Qii 


Qn 


Qi 2 

Q22 







Qn 

Q21 

Qli 


Qli 


QH 



Figure 1: Macro elements: = Qj^k U Qj,k+l U Qj+i,k U Qj+i,fc+i, (J, K) = (j, k) 

We will employ capital letters to indicate odd integer indices for those macro 
patterns on the macro element. Owing to Theorem 2.1, the global checker-board 
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pattern basis function Ch in (2.4) can be expressed explicitly as follows: 

c. = (2.5) 

JK 

We now try to stabilize x minimally so that the modified pairs 

fulfill the uniform inf-sup condition. In this section we introduce the stabiliza¬ 
tion of pressure approximation space by eliminating one-dimensional global 
checker-board patterns from . Alternatively, the stabilization of velocity ap¬ 
proximation space again with a globally one-dimensional modification, 

is given in §3. 

2.3. Stabilization of 

Define as the L^(f2)-orthogonal complement of in £Pq, that is, 

^1} = © ^, dim (^) =Nq-2. (2.6) 

We are now ready to propose our Stokes element pair as follows: 

^ ‘i™(^"o'‘ x^o) = + Nq-2. (2.7) 

2 . 4 . The discrete Stokes problem 

Now define the discrete weak formulation of (2.2) to find a pair {xih^Ph) S 
X such that 

ah{uh,^h) - bh{-Vh,Ph) = (f,v/i) (2-8a) 

bh{uh,qh) = 0 (2.8b) 

where the discrete bilinear forms ah{-,-) ■ ^"0^ x K and bh{-,-) : 

X —>• ffi. are defined in the standard fashion: 

Nq Nq 

a/i(u,v) = z/^(Vu, Vv)q^ and bh{v,q) ■ v,q)Q.. 

j=i i=i 

As usual, let | • |i,/t denote the (broken) energy semi-norm given by 

|v|i,/» = Va/i(v,v), 
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which is equivalent to || ■ on ^”g^. Also, denote by || • \\m,h and | • \m,h the 
usual mesh-dependent norm and semi-norm; 




1/2 

■ 

l|v|U.h = 

ll''llff-(Q) 

and |v|m,;i = 



■ 


■ QG^k 


respectively. Let denote the divergence-free subspace of to i.e., 

= {v/j e I bh{vh,qh) = 0, yqh & ^q}- (2.9) 

Then the solution u^, of (2.8) lies in and satisfies 

ah{uh,Vh) = (f,v/,) Vv/, e (2.10) 

We state the main theorem of the paper, whose proof will be given in §4. 
Theorem 2.2. g^ X satisfies the uniform discrete inf-sup condition: 

sup > /3||g^||g ^ Vg?, € (2.11) 


3. Alternative stabilization by enriching the velocity space 

In this section we consider an enrichment of by adding a global one¬ 

dimensional bubble function space based on the quadrilateral nonconforming 
bubble function [5, 6, 8, 15]. We then compare two proposed nonconforming 

finite element space pairs ^"g^ x .^g and ^" g^ x .^g. Indeed, these two 
spaces very closely related. The velocity solutions obtained by these two spaces 
are identical while the difference between the two pressures isof order 

On a reference domain Q := [—1,1]^, the DSSY nonconforming element 
space is defined by 

DSSY{Q) = Span{l,a;,y, 0fc(x) -Okiv)}, 


where 


dk{t) 


J — |t^, fc = 1, 
\t^ -ft^ + lt^, k = 2. 
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Let Fq ■. Q ^ Q he d. bijective af&ne transformation from the reference domain 
onto a rectangle Q. Then define 


DSSY{Q) = {r; o v€ DSSY{Q)^ . 


(3.1) 


The main characteristic of DSSY{Q) is the edge-mean-value property: 

Ip ds = ■j/>(midpoint of E) Vi/) € DSSY{Q), (3.2) 

where denotes [8, 15]. 

The vector-valued DSSY nonconforming finite element space is defined by 

DSSYj = {v € I V, := vjg^ € DSSY(g,) Vj = 1, • • • ,Nq-, 

V is continuous at the midpoint of each interior edge 

and vanishes at the midpoint of each boundary edge in ,^h}. 

For each macro element define £ DSSYq such that supp(i/)Q^) C 

and its integral averages over the edges in vanish except on the two 
edges dQj^i fl dQj+i^t, £ = k,k + 1 : 


f ds = u, f ipQ 

d QQj,ke\dQj^\^k d dQj,k+ie\dQj + i,k+i 


M ds = —V. 

JK 


where u denotes the unit outward normal vector of Qj^i on the edge dQj^i fl 
cig^+i^^, £ = k,k F Define a basis function for the global bubble function, as 
shown in Figure 2, and a space of global bubble functions as follows: 


= Span {B^i}, 




(3.3) 




We are now ready to enrich as follows: 


^nc,h _ ^nc,h 

^ 1,0 ~ *^ 1,0 


(3.4) 


Remark 3.1. The dimension of the pair of spaces .^”g x is 2iV® -|- Nq. 

We state the uniform inf-sup stability as in the following theorem, whose 
proof will be given in §4. 


Theorem 3.2. g^ X tY’l} satisfies the uniform discrete inf-sup condition. 
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"^■+i,fc+i 

Qj.,k+1 

Qj+1,A;+1 


' - O 




brj./c 

Qj,k 

Qj+i,fc 


f 


1 ™7+J,ic-l 


Figure 2: The basis function '^PqM € DSSYq, associated with the macro element 

takes the value u and —u at the midpoints m. . _ i and m. . , i , respectively, and value 0 at 

2 2 

any other midpoints m’s shown in the figure. = Qj,k hJ Qj,k+i U Qj+i,fc U 


3.1. Comparison between x 3^^ and x 

In this subsection, we will compare the two nonconforming finite element 
space pairs x .^g and .^"g^ x .^^g . These two pairs are closely related 

such that X ^g can be understood as a slight modification of g^ x ^g . 

For ,^"g^ X i3^g , we have the following discrete weak formulation: Find a 


pair {u'i^,p'i^) G .^"g^ x ^g such that 

ah{K,Vh)-bh{v',^,p'h) = Vv;^ € (3.5a) 

bh{K,qh) = 0 e ^0 . (3.5b) 

Let denote the divergence-free subspace of .^"g^ to i.e., 

= {v; e I 6.(v', q’^) = 0, Vg; e J^gh- (3.6) 

Then the solution of (3.5) lies in and satisfies 

a/i(ub<) = (f,v;) Vv^ e (3.7) 
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The following lemma implies that the two divergence-free subspaces defined in 
(2.9) and (3.6) are identical, that is, our two proposed nonconforming finite 

element space pairs x and x produce an identical solution 

for velocity. 

Lemma 3.3. The spaces and defined by (2.9) and (3.6), respectively, 

are equal. 

Proof. Let v/i € be given. Since g), G Span{.^g 0'^^} and by Theorem 2.1, 
we get bh{vh, q^) = 0. This implies Vh G so C It remains to prove 
Q i^g given, where v/h G and G 

In particular, if we consider G then bhi'^'h, q'ff) — 0 implies b;i = 0. 
Therefore G and bh{'v'f^,qh) = 0 for any qu G 3^^ since 3^^ C .^g- 

Hence G which shows C . This completes the proof. □ 

Owing to Lemma 3.3, where U/j and w.'^ are the solutions of (2.8) 

and (3.5), respectively. Moreover, the difference between the two pressure solu¬ 
tions obtained by (2.8a) and (3.5a) fulfills 

bhi-Vh^p'h -Ph) =0, Vv/i G ^1,0^. 

By Theorem 2.1, — Ph & that is, p), can be represented by 

Ph =Ph + a<Ch, a G K. 

Taking =Bh & 3§^ in (3.5a), we obtain 

ahh{Bh,Ch) = ah{uhi^h) — - bhiBh,Ph), 

Nq 

= Uh,'VMh)Q^ - (f,B/t) -bhiBh,Ph), 

1=1 

/du \ 

= l^y^(-Au/t,B^)Q. +v( -^,Bh ) - (f,B/j) - 6;j(B/j)|^;^) 

j=l \ ! dQi 

Since the solution is a piecewise linear polynomial, that is, G ^"g^, the 
first term in (3.8) is equal to zero. And we easily check that the second and 
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last terms in (3.8) turn out to vanish by the characteristics of the space . A 
simple calculus using the Divergence Theorem yields 

C/j) =-. (3.9) 


Invoking (3.9), one obtains 

(f,B/.) 

bh{^h,'Ch) 

Hence, p'^-ph = -h{f,Mh)Ch. 




(3.10) 


We summarize the above result as follows: 

Theorem 3.4. Let {uh,Ph) G x and (u'^,p),) G x are the 

solutions of (2.8) and (3.5), respectively. Then 


= u), and ph - p'h = -h{f,'Bh)Ch. 


(3.11) 


3.2. Interpolation operator and conference results 

We recall from [21] that the global interpolation operator XI/i : H^(H) — 

^nc,h 

is defined through the local interpolation operator IIq : H^(Q) —>■ 

^nc.?i(Q) 

n;,|Q = nQ yQetTh. 

Here, IIq is explicitly defined by 

nQw(Mfe) = + ^ (3^^2) 

where 14 -1 and 14 are the two vertices of the edge Ek with midpoint Mk of Q. 
Define an interpolation operator St '■ H^{Lt) n To(^) ^ t>y 

{Shq,z) = {q,z) Vze^l}. 

Since 11?, and Sh reproduce linear and constant functions on each element 
Qj G and macro element respectively, the standard polynomial ap¬ 

proximation results imply that 

jjv - n?iv||o -k hjv - + h'^\v - n?,v| 2 ,f, (3.13a) 

+/ii/2|v - n?,v|o.aa < Ch^hh Vv e H2(0), 
|| 9 -^h<?||o,n <C/i||g||i VqGH\n)nLl{n). (3.13b) 
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Owing to (3.13), a standard application of Theorems 2.2 and 3.2, and the 
second Strang lemma yields the following optimal error estimate: 

Theorem 3.5. Assume that (2.1) is H‘^{D,)-regular. Let (u,p) and {uh,Ph) 
be the solutions of (2.2) and (2.8) (or (3.5)^ respectively. Then the following 
optimal-order error estimate holds: 

||u - Uh\\o + h [|u - Uh\i,h + \\p-Ph\\o] < Ch"^{\u \2 + Iblli). 

Remark 3.6. In the above theorem, after the result for x is shown, 

the corresponding result for x to Theorem 3.5 can be obtained a 

combination of Theorems 3.4 and 3.5. The order of two spaces x 

and g^ X can be of course exchanged. 


4. Proofs of Theorems 2.2 and 3.2 

In this section we will show that .^"g^ x J^g and .^”g^ x 3^^ satisfy the 
uniform discrete inf-sup condition. For this, some useful results [10, 23] will be 
used; in particular. Lemma 4.1, a result of Qin [23], will be utilized. 

Our proof starts with setting 



Then denote by the L^(fl)-orthogonal complement of in such that 


^’4, dim {^(f) =Nq-2 and dim {W^) = ^Nq - 1. (4.1) 
Let denote the discrete divergence-free subspace of .^"g^ to 3^(f, that 


is. 


Z'* 



bhiyh^qh) = 0\/qh G } ■ 


Considering the conforming bilinear element 


= {v,h€ Hj(L!) I each component of 'v^hlg is bilinear VQ G , (4.2) 
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and denote the discrete divergence-free subspace of .^i q to that is 
= {v,„ e I bh{vi,h,qh) = 0 V% e 3^^} . 


Denote by S‘2h and ^ 2 ?i the sets of all edges and interior edges, respectively, in 
. Set to be the subspace of defined by 


^me,2h 

*^1.0 


V/i 


G ^ 


nc,h 

1,0 




r^e-?2fe \or"/ 




(4.3) 


where ipr^ € is the basis function associated with the midpoint of the 

macro edge G as described in detail in the caption of Figure 3. Notice 
that = Ni-l. 



Figure 3: The basis function Tp^M € q^, associated with the macro edge F^ = X’ 

takes value 1 along the four line segments joining the midpoints rrijj^s fc+i ’ k ’ 

and fe—l; s,nd value 0 at any other midpoints m’s shown in the figure. denotes 

the midpoint of the macro edge F^^^ the common edge of the two macro elements 
and Qj_|_ 2 ,x’ Qj,k ~ Qj,k-\-i U Qj-\-i,k U S’Hd Qj^ 2 ,k ~ Qj+ 2 ,k U 

Qj^2,k+l U Qj+3,fc U Qj+3,fc+l- 


Next, we quote the Subspace Theorem of Qin as in the following lemma: 

Lemma 4.1 ([23]). Given V^ x , let Vi and V 2 be two subspaces ofY^ and 
Pi and P2 be two subspaces of P^. Let the following four conditions hold: 
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(1) ph =Pi+P2; 

(2) there exist j3j > 0,j = 1,2, independent of h, such that 



(3) there exist Uj > 0, j = 1, 2, such that 

\bhivj,qk)\ < ctjWjli.hhkWop, \/Vj eVj and Vqk e Pk, j,k = 1,2] j k, 


with 


aia2 < /3 i/32. 


Then, x satisfies the inf-sup condition with the inf-sup constant depending 


only on ai,a 2 , Pi, ^ 2 - 
4-1. Proof of Theorem 2.2 

The following lemma is an immediate consequence of the Divergence Theo¬ 
rem, which will be useful to prove Lemma 4.3: 

Lemma 4.2. Let Q C be a rectangular domain. Suppose that w is a two- 
variable function whose components are bilinear polynomials on Q. Then the 
following holds: 



Lemma 4.3. x satisfies the uniform discrete inf-sup condition: 



(4.4) 


Proof. We begin with invoking [4] that x satisfies the uniform inf-sup 
condition, that is, there exists a positive constant P independent of h such that 



(4.5) 
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Let Qh G W^,qh 0 be arbitrary. Then, (4.5) is equivalent (cf. [10], p. 118) to 
the existence of Vf,?i S such that 

bh{'Vbh:qh) = \\qh\\l.n: (4-6a) 

|vbh|i,n < ^ll^llo.n- (4.6b) 

Now Lemma 4.2 implies that Ilh'Vbh G Z^ and 


bh(nhVbh,qh) = bh{vbh,qh) = lk?i|lo.n- 


(4.7) 


By Young’s inequality, the definition of interpolation operator 11^ and (4.6b), 
one sees that 

< C|vb?i|i_n < -^\\qh\\o,Q, (4.8) 

where the constant C is independent of mesh size h. Notice that the element 
of v/i = Il/iVb/i G Z^ satisfying (4.7) and (4.8) plays a role of an equivalent 
statement to (4.4). Hence the lemma is complete. □ 


Lemma 4.4. x satisfies the uniform discrete inf-sup condition: 

bh{'Vh,qh) 


sup 

/TATne ,2h 


\^h\l,h 


>^lk/i||o,a 'iqh&^c 


(4.9) 


Proof. Set 

= {q€ Llin) I q\QM e PoiQ^) G dim (J^^^) = Nq/A - 1. 


Due to Lemma 3.1 in [22], x satisfies the uniform inf-sup condition, 

that is, there exists a positive constant /3 independent of h such that 


Let qh 

where 


sup Vg,G.^o^". (4.10) 

= ctJKp’jK € be arbitrary. Consider = ''^ajxPJK G 

JK JK 

Pjk = Xqy^- Then there exists v/i = ( ^ ) ^ 

\br^ J 
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such that (4.10) holds. From this Vh, we define Vh € ^ as follows: 


v^ = - E 

r"G<f2h 



Then the following three equalities are obvious: 


||9?i||o,n = ll9/illo,a, 

bh{'Vh,qh) = bh{vh,qh). 


(4.11a) 

(4.11b) 

(4.11c) 


From (4.10) and (4.11), the inf-sup condition (4.9) for x follows. 

This proves our assertion. □ 


Utilizing Lemma 4.1, we are ready to prove Theorem 2.2. 


Proof of Theorem 2.2. We will check the conditions of Lemma 4.1. Let Vi = 
Z^,V 2 = Obviously, Y, and Pp j = 1,2 are 

subspaces of and j respectively, so that Condition (1) holds. Moreover, 

Lemmas 4.3 and 4.4 imply that Condition (2) holds. Since bh{'Vi,q 2 ) = 0 holds 
for any vi G Vi and any q 2 G P 2 , one has ai = 0. Consequently, Condition 
(3) holds. Hence by Lemma 4.1, x satisfies the inf-sup condition 

( 2 . 11 ). □ 


4 . 2 . Proof of Theorem 3.2 

In order to prove Theorem 3.2, the following lemma is needed. 

Lemma 4.5. x satisfies the inf-sup condition, that is, there exists a 
positive constant /3 independent of h such that 

sup > y 3 ||qr^||g (4.12) 

Proof. Let qh G be given by qh = aCu with a constant a G M, and set 
v/i = B/i G Recall (3.9) so that 

bhivh,qh) = Y- (4.13) 

a 
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Also, it is trivial to see 


ll9/i||o,n = |a|. 


(4.14) 


It remains to compute For this, we notice that |v?t|i_Q does not depend 

on the mesh size h of Q, since it is a two dimensional region. Indeed, there exists 

_^ r (j 

a constant Ci independent of h such that |v/i|f ^ = ^ L |V V/Ip dx = 

ITT J 0. 

Hence, we get 




c 


Wh\i,h = where C = 

h 


(4.15) 


Now, the combination of (4.13), (4.14) and (4.15) leads to (4.12) with the inf-sup 
constant /3 = 1/C. This completes the proof. □ 


Theorem 3.2 is now ready to be shown, by using Lemmas 4.1 and 4.5. 

Proof of Theorem 3.2. Let Vi = V 2 = and = ;^, P 2 = 

Since 0 Condition (1) in Lemma 4.1 holds. Moreover, The¬ 

orem 2.2 and Lemma 4.5 imply Condition (2) holds. Finally, bh{'vi,q 2 ) = 0 
holds for any vi S Vi and q 2 € P 2 by Theorem 2.1, which implies that ai = 0. 

Consequently, Condition (3) holds. Hence, <^"g^ x satisfies the inf-sup 
condition. Note that the constant in each step is independent of h. □ 


5. Numerical results 

Now we illustrate a numerical example for the stationary Stokes problem on 
uniform meshes on the domain ft = (0,1)^. Throughout this numerical study, 
we fix = 1. 

First we calculate the discrete inf-sup constants of various finite element 
pairs including our suggestions. 

In contrast to the ^(/i)-dependent inf-sup constant of conforming bilinear 
and piecewise constant finite element pair [3, 4], our two proposed noncon¬ 
forming finite elements satisfy the uniform inf-sup condition at least on square 
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h 

Pi 

Order 

P2 

Order 

Ps 

Order 

1/4 

4.9642E-01 

- 

4.9560E-01 

- 

5.0000E-01 

- 

1/8 

2.8605E-01 

0.78 

4.6791E-01 

0.08 

4.6746E-01 

0.09 

1/16 

1.5029E-01 

0.93 

4.4415E-01 

0.07 

4.5296E-01 

0.04 

1/32 

7.6544E-02 

0.97 

4.2863E-01 

0.05 

4.4526E-01 

0.02 

1/64 

3.8562E-02 

0.99 

4.1864E-01 

0.03 

4.4051E-01 

0.02 


Table 1: Estimation of the inf-sup constants = 1,2,3, for the three finite element pairs 

V £anc,h ^ric,h 

•^ 1,0 ^ ^0 ’ ^ 1,0 ^0 ’ ^ 1,0 ^0 ■ 

meshes. To confirm theoretical analysis, we give the numerical results of the 
discrete inf-sup constants [19] in Table 1. 

We will borrow the two numerical examples from [22] . The source term f is 
generated by the choice of the exact solution. 

u(x, y) = {s{x)s'{y), -s(y)s'(x)), p(x, y) = sm{2nx)f{y), (5.1) 

where s{t) = sin(27rt)(t^ — t) and s'{t) denotes its derivative. The velocity u 
vanishes on dQ and the pressure p has mean value zero regardless of /. 

Several interesting numerical results for the pair x are presented, 

while the corresponding numerical results for the pair x are omitted 

here, since they behave quite similarly to those case for the pair x . 

Numerical results with f{y) = 3 _^,hii y are shown in Table 2. We observe opti¬ 
mal order of convergence in both velocity and pressure variables. Also numerical 
experiments are carried out and presented in (5.1) for f[y) = 25 -io\an^ y 
which has a huge slope near the boundary on y = 1. Since the pressure changes 
rapidly on the boundary y = 1, convergence rates show a poor approximation 
in coarse meshes in Table 3. However, as the meshes get finer, optimal order 
convergence is observed as expected from the inf-sup condition. 

The following numerical results highlight the reliability of our proposed finite 
element space compared to the case of using the conforming bilinear element for 
the approximation of the velocity field. Recall that the pair of conforming finite 
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h 

|u - 

Order 

l|u- wllo 

Order 

\\p-Phh 

Order 

1/4 

1.5087E-0 

- 

2.1583E-1 

- 

2.2190E-1 

- 

1/8 

8.1269E-1 

0.8926 

5.5033E-2 

1.9715 

1.4098E-1 

0.6544 

1/16 

4.1360E-1 

0.9745 

1.3930E-2 

1.9821 

6.4738E-2 

1.1229 

1/32 

2.0767E-1 

0.9939 

3.4936E-3 

1.9954 

3.2509E-2 

0.9938 

1/64 

1.0394E-1 

0.9985 

8.7411E-4 

1.9988 

1.6411E-2 

0.9862 

1/128 

5.1985E-2 

0.9996 

2.1857E-4 

1.9997 

8.2359E-3 

0.9947 

1/256 

2.5994E-2 

0.9999 

5.4646E-5 

1.9999 

4.1222E-3 

0.9985 

1/512 

1.2997E-2 

1.0000 

1.3661E-5 

2.0000 

2.0616E-3 

0.9996 

1/1024 

6.4987E-3 

1.0000 

3.4154E-6 

2.0000 

1.0309E-3 

0.9999 


Table 2: Numerical results for uniform meshes with f(y) = - —-— 2 — 

o—XEiii y 


element space combined with the piecewise constant element space x 
is unstable unless f is smooth enough as quoted in the following Corollary: 

Corollary 5.1 (Boland and Nicolaides, Cor. 6.1 in [4]). For /3 G (0,1), there 
exists f G L^(f2) such that the pressure approximation to (2.2) by using ^J’g x 
fulfills 

lb-P/i||o > C'h'^||f||o forh<hp (5.2) 

for some hp > 0, independent of h. 

With P = 0.3 fixed, some comparative numerical results for conforming and 
nonconforming pairs using .S^’g x and ^”g^ x are shown in Tables 4 
and 5, respectively. These results ensure the superiority of our nonconforming 
method over the conforming counterpart. 

Throughout our numerical experiments, the 4x4 Gauss quadrature rule is 
adopted for each rectangular element. The approximate data for f are calculated 
by following the proof of Theorem 6.1 in [4] at the 4x4 Gauss points in each 
element of 512 x 512 mesh. The reference solutions used in error calculation are 
obtained by using the DSSY element [8] with the 512 x 512 mesh. The graphs 
of components of f are given in Figure 4. 
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h 

|u - 

Order 

l|u- U/^||o 

Order 

lb-Philo 

Order 

1/4 

1.5086E-0 

- 

2.1578E-1 

- 

1.7459E-1 

- 

1/8 

8.1268E-1 

0.8925 

5.5016E-2 

1.9716 

1.1835E-1 

0.5609 

1/16 

4.1360E-1 

0.9744 

1.3926E-2 

1.9820 

5.7158E-2 

1.0501 

1/32 

2.0767E-1 

0.9939 

3.4938E-3 

1.9950 

3.6347E-2 

0.6531 

1/64 

1.0394E-1 

0.9985 

8.7450E-4 

1.9983 

2.3178E-2 

0.6491 

1/128 

5.1985E-2 

0.9996 

2.1872E-4 

1.9993 

1.3569E-2 

0.7725 

1/256 

2.5994E-2 

0.9999 

5.4690E-5 

1.9998 

7.3091E-3 

0.8925 

1/512 

1.2997E-2 

1.0000 

1.3673E-5 

1.9999 

3.7516E-3 

0.9622 

1/1024 

6.4987E-3 

1.0000 

3.4183E-6 

2.0000 

1.8899E-3 

0.9892 


Table 3: Numerical results for uniform meshes with f(y) = —— n-h 

25 —lOtan-^y lU 

Remark 5.2. It should be stressed that the degrees of freedom for both 
and X are essentially identical, although numerical results are guite 

different. Further investigations need to be sought to analyze the differences 
between the conforming bilinear element and the Pi nonconforming element. 


The X component o11(x,y) 

X10 ; 



The y component o1 f(x,y) 



Figure 4: The graph of data f 
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h 

l^re/ 

order 

llUre/ -U/ill 

0 order 

llPre/ Philo 

order 

1/4 

2.8248E-2 

- 

1.8470E-3 

- 

7.2967E-2 

- 

1/8 

1.6008E-2 

0.8193 

5.3114E-4 

1.7981 

5.6105E-2 

0.3791 

1/16 

8.5909E-3 

0.8980 

1.4266E-4 

1.8964 

4.1920E-2 

0.4205 

1/32 

4.4824E-3 

0.9385 

3.7531E-5 

1.9265 

3.1925E-2 

0.3929 

1/64 

2.3084E-3 

0.9573 

9.6932E-6 

1.9531 

2.4932E-2 

0.3567 

1/128 

1.1939E-3 

0.9512 

2.4703E-6 

1.9722 

1.9829E-2 

0.3304 

1/256 

6.4542E-4 

0.8874 

6.2940E-7 

1.9727 

1.5938E-2 

0.3152 


Table 4: Numerical results for X 

^g when (9 

= 0.3 


h 


order 

||u^e/ - U,i|| 

0 order 

||Pre/ ~ Ph ||o 

order 

1/4 

2.8359E-2 

- 

1.8561E-3 

- 

4.9406E-2 

- 

1/8 

1.7966E-2 

0.6585 

5.0224E-4 

1.8858 

2.6963E-2 

0.8737 

1/16 

1.0379E-2 

0.7916 

1.3390E-4 

1.9072 

1.4305E-2 

0.9144 

1/32 

5.6226E-3 

0.8844 

3.5144E-5 

1.9298 

7.5726E-3 

0.9177 

1/64 

2.9406E-3 

0.9351 

9.0617E-6 

1.9554 

3.9235E-3 

0.9486 

1/128 

1.5002E-3 

0.9710 

2.3029E-6 

1.9763 

1.9663E-3 

0.9966 

1/256 

7.3601E-4 

1.0274 

5.7096E-7 

2.0120 

8.9372E-4 

1.1376 


Table 5: Numerical results for x when /9 = 0.3 
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